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Abstract 

We study the transport properties of a superconducting Nb film with a square lattice of artificial 
pinning centers (APCs) as a function of dc current, at a temperature close to the superconducting 
transition temperature of the film. We find that, at low dc currents, the differential resistance of 
the film shows the standard matching field anomaly, that is, the differential resistance has a local 
minimum at magnetic fields corresponding to an integer number of flux lines per APC. However, 
at higher dc currents, the differential resistance at each matching field turns to a local maximum, 
which is exactly opposite to the low current behavior. This novel effect might indicate that the 
flux lines in the APC system change their flow mode as the dc current is increased. 
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Superconducting films with a lattice of artificial pinning centers (APCs) have been in- 
vestigated for several yeara^ i^i^i^i^ . It is found that the vortex lattice is highly ordered at 
certain matching fields, where the number of flux lines equals an integer multiple n 
of the number of APCs. This ordering effect provides additional pinning interactions that 
result in a local minimum resistance and maximum critical current at the matching fields. 
Interest has been growing in such systems in the presence of driving forces^ . Although much 
work has been done in systems with random disorder^^, the behavior of AFC systems is 
only beginning to be explored. In this letter, we investigate the transport properties of a 
superconducting Nb film with a square lattice of APCs as a function of dc current through 
the film. Our results show novel behavior in a specific range of the dc current, which has 
not been observed in previous experiments. 

Figure 1(a) shows a schematic of the device discussed in this letter. This device was 
patterned onto an oxidized Si substrate using laser interferometric lithography. The entire 
device is made of a single layer of Nb film (100 nm thick) with a square lattice of APCs with 
period d^l fim and diameter fim. A scanning electron microscopy (SEM) image of 

the Nb film is shown in Fig. 1(b) with a 2 /xm size bar. Technical details of the patterning 
procedure can be found in Ref. [3]. 

Figure 2 shows the differential resistance of the Nb film as a function of magnetic field 
normal to the film with different dc bias currents (J^c) through the sample. The resistance 
was measured by conventional four-terminal techniques using an ac resistance bridge. The 
temperature of the system was fixed at T=7.07 K, close to the superconducting transition 
temperature, Tc=7.13 K, of the film. All the resistance curves show oscillations as a function 
of magnetic field, with a fundamental period of ~ 18 G, corresponding to one superconduct- 
ing flux quantum h/2e per unit cell of the square lattice, or alternately, one flux line per 
APC. At low dc bias currents (/dc<100 /iA), the differential resistance curve has a local 
minimum at magnetic fields Hn corresponding to an integer number n of flux lines per APC. 
This behavior is the one that is conventionally observed. It is associated with strong pinning 
of flux lines at the matching fields Hn, leading to lower dissipation associated with flux line 
motion, and hence a lower resistance. At higher dc bias currents (/dc^HO /iA), however, the 
local minimum becomes a maximum at each matching field, exactly opposite to the results 
at the low dc current. The current at which this change occurs is between 100 fiA and 110 
/iA, where the behavior of the differential magnetoresistance shows a dramatic change in a 
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small current range of 10 /iA. At even higher dc bias currents (/dc>160 /iA), not plotted out 
in Fig. 2, the magnetoresistance oscillations are washed out. 

The novel behavior observed in Fig. 2 can also be seen in Fig. 3(a), the differential 
resistance measurements as a function of dc current at two different magnetic fields of 18 
G and 9 G, corresponding to the first matching field Hi, and Hi/2, respectively. This 
measurement has been done by the same four-terminal technique used for Fig. 2, and at 
the same temperature T=7.07 K. Figure 3 illustrates that the dV/dl vs. Idc behavior can 
be divided into three regimes: I, II and III^. In regime I, the measurements show zero 
differential resistance, which indicates that the vortex lattice is static due to the strong 
pinning force. In regime II, the differential resistance increases gradually when we increase 
the dc current. This regime is called the plastic flow regime, where the vortex lattice has been 
partly disordered and shows plastic flow behavior. In regime III, the measurements show a 
sharp resistance jump to the normal state of the film, which indicates that the vortex lattice 
has been completely disordered and shows avalanche motion. In our measurements, there 
exist in addition three sub-regimes in the plastic flow regime, which are labeled 11(1), 11(2) 
and 11(3) in Fig. 3. In 11(1) the measurements show the standard matching field anomaly, 
i.e., the 18 G curve, where each flux line is expected to be tied to one APC, has a smaller 
resistance compared with the 9 G curve. In 11(2), the differential resistance measured at 18 
G is larger than the one measured at 9 G, which is opposite to the results in regime 11(1). 
In 11(3), the 18 G curve and the 9 G curve merge together. No oscillations appear in this 
sub-regime. It should also be pointed out that the crossing point between regime 11(1) and 
11(2) is quite sensitive to the system temperature. As shown in Fig. 3(b), this crossing point 
moves from 103 /iA to 152 /iA when we cool our system down just 0.01 K, from 7.07 K to 
7.06 K. In other words, the novel behavior in the regime 11(2) only happens in a very small 
temperature range, which might be the reason why it has not been observed in previous 
work. 

To our knowledge, no current theory can explain this unusual behavior. We propose 
here a phenomenological picture to explain our experimental results. The total force / 
acting on each flux line in a superconducting film with APCs can be expressed as / = 
jfp _j_ _|_ _|_ f ^ where f"^ is the pinning force due to the trapping of the flux hne 
by a APC, /™ is the repulsive interaction with other flux lines, frp the force associated 
with thermal fluctuations, and the Lorentz force, present when transport current is sent 
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though the film^. The first two terms in / would result in a highly ordered vortex lattice 
configuration, especially at each matching field; the last two terms drive the system from 
the highly ordered vortex lattice state to a disordered state. The motion of fiux lines is 
determined by the balance between f"^ + f^^ and frp + f^. In our experiments, we fix the 
system temperature at a value slightly below the superconducting transition temperature of 
the Nb film, that is, we set fj^ to be constant, and vary by tuning the dc current through 
the device. 

The behavior demonstrated in Fig. 3 might be associated with a change of the fiow mode 
of the fiux lines when the dc current is increased. At very low dc currents (regime I), the 
pinning force dominates the system, {f^^ + f^^)^{fT + fd)- -^^ ^ consequence, the fiux lines 
are static and the measurements show zero differential resistance. This corresponds to the 
situation shown in Fig. 4(a), which shows a schematic of the potential landscape seen by 
the vortex lines at Hi. There is a minimum of potential at each APC that pins exactly one 
vortex line. Any vortex line that is displaced from its pinning center by thermal fluctuations 
is immediately forced back due to the repulsive interactions between vortices, leading to a 
minimum in the resistance due to vortex motion. For half-fllling (corresponding to Hi/2, 
the situation shown in Fig. 4(b)), the distance between the vortices is larger, leading to a 
smaller vortex-vortex interaction. Each vortex is therefore not so rigidly pinned as in Fig. 
4(a), giving rise to a larger resistance due to vortex motion. 

The effect of applying a dc current is to add a flnite slope to the potential proflle, as shown 
in Figs. 4(c)-(f). The larger the current, the steeper the slope. At low I^c, (corresponding 
to region 11(1)), the pinning potential is still strong enough that the situation described 
above for Figs. 4(a) and (b) still holds, and the resistance at Hi (corresponding to Fig. 
4(c)) is less than the resistance at Hi/2 (Fig. 4(d)). At still higher I^c (corresponding to 
region 11(2)), however, the situation changes. Consider the case of full fllling, at Hi, shown 
schematically in Fig. 4(e). Here, the effective pinning potential at one end of a pinning 
center is effectively reduced by the slope of the potential induced by the dc fleld, so that 
thermal fluctuations can easily kick a flux line out of the potential well. If this happens in 
the center of the superconducting fllm, the repulsive interactions between vortices are still 
likely to push the vortex line back into its pinning site. If the vortex pops out of a pinning 
site at the edge of the fllm, however, it will leave the fllm altogether. At this point, the 
flux line at the pinning site immediately adjacent to the now empty site will no longer feel 
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the repulsive interaction due to the absent vortex hne. If the repulsive interaction from its 
remaining vortex neighbors is sufficient to overcome the reduced pinning potential, it will 
jump to the empty pinning center. This jump mechanism will be continued down the linear 
chain of vortices, leading to a mass movement of vortices in the lattice, and hence a higher 
resistance. For half-filling {Hi/2, shown in Fig. 4(f)), the smaller repulsive interaction 
between vortices, arising from the greater inter-vortex distance, is not sufficient to overcome 
the effective pinning potential, so that the flux lines do not move. Hence the resistance 
of the sample at Hi/2 is lower than its resistance at Hi. This mechanism clearly depends 
on a flne balance between the dynamic pinning potential, which depends on Idc, and the 
repulsive vortex-vortex interactions, and hence can only be seen in a narrow range of Idc- 
At even higher currents, corresponding to region 11(3), the potential proflle is tilted so much 
by the applied dc current that the vortex lattice flows at both Hi/2 and Hi, so that no 
clear oscillations are seen in the magnetoresistance. A similar picture also holds at higher 
magnetic fleld, where there is more than one vortex per unit cell. One might also point out 
that vortex-antivortex pairs are present even at zero fleld, as evidenced by the presence of 
a flnite resistance. 

An alternative explanation is based on the remarkable similarity between the data shown 
in Fig. 2, and the oscillations in conductance observed in an Aharonov-Bohm interferometer 
formed from a quantum dot^^. In that experiment, a similar change in the phase of the 
oscillations is observed as the phase difference between the two arms of the interferometer is 
varied by means of an external gate potential. Indeed, in that experiment, the oscillations 
at the fundamental period $o disappear at the cross-over point, and only oscillations at 
$0/2 remain. Hints of this behavior are also seen in our data (see curve (d) in Fig. 2). This 
suggests that the transport current in our system is modulating the phase of the interference 
around each antidot in the lattice. The exact nature of this modulation is not clear to us at 
the moment, and further investigation is required to understand the origin of this unusual 
behavior. 

This work is supported by the NSF under grant number DMR-0201530. 
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Fig. 1, Z. Jiang et al. 




FIG. 1: (a) Schematic of the device. The entire device is made of a single layer of Nb film with a 
square lattice of APCs. The middle section between the two voltage leads, the part measured, is 
50-/xm-wide and 290-/xm-long. (b) A scanning electron microscopy (SEM) image of the Nb film. 
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Fig. 2, Z. Jiang et al. 
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FIG. 2: Differential resistance of the Nb film as a function of magnetic field normal to the film 
at different dc bias currents (/^c) a temperature of r=7.07 K: (a)/cic=61 yuA; (b)J(ic=70 /xA; 
(c)7de=100 /xA; (d)Jde=110 M; (e)/dc=120 M; (f)^dc=140 M; (g)^dc=159 M- 
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Fig. 3, Z. Jiang et al. 
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FIG. 3: Differential resistance of the Nb film as a function of dc current at two different magnetic 
fields of 18 G (solid line) and 9 G (dashed line), at two different temperatures: (a) r=7.07 K and 
(b) r=7.06 K. 
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Fig. 4, Z. Jiang et al. 




FIG. 4: Schematics of the potential landscape seen by the vortex lines at different dc bias currents 
and at two different magnetic fields: Hi (left column panel) and Hi/ 2 (right column panel). Each 
flux line is denoted by an arrow. As predicted by theory, an empty AFC shows a deeper potential 
compared with one filled by a flux lineiS. 
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